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^ . Abstract 

Phase structure and phase transitions in dense QCD are studied using the 
Schwinger-Dyson (SD) method in the improved ladder approximation. We construct 
the Cornwall- Jackiw-Tomboulis (CJT) effective potentials at finite temperature for two 
types of pairing ansatz, namely the Color-Flavor locking (CFL) state and the two flavor 



superconducting (2SC) state. Strong coupling effects at low densities, such as the off- 
Fermi surface and antiquarks contribution to the pairing gap due to the large effective 
coupling, make the gap, critical temperature and those ratio deviate from those weak 
coupling values. Nevertheless, the ratio of the physical quantity for the CFL and that 
for the 2SC does not differ so much from the weak coupling value, as long as the pairing 
interaction is taken to be the same for the both cases. As a consequence, the CFL state 
always dominates over the 2SC state and the critical temperatures to the quark-gluon 
plasma (QGP) phase from both states coincides in the chiral limit. The energy gain in 
the CFL state relative to the 2SC state gets smaller towards the critical line dividing 
(/i, T) plane into the QGP and CFL phases, and thus small perturbations can remove 
the degeneracy of these critical lines. As one of such perturbations, the effect of the 
strange quark mass m s on the quark-pairing is examined. In particular, using a simple 
kinematical criterion, we discuss the behaviour of the " color-flavor unlocking line" on 
which the CFL phase turns into the 2SC phase, against the variation of m s . 
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1. Introduction 



It is believed that the QCD ground state changes against the variation of the external 
conditions, such as the baryon density, the isospin density and the strange quark mass. In 
particular, the color superconducting phases at high baryon density have attracted much 
interest in high density QCD. The asymptotic freedom of the QCD ^ 2 ) leads to the smaller 
effective coupling at higher densities : This gives us the conjecture of deconfined quark 
matter 3 ) in which quarks and gluons constitute the active degrees of freedom and are weakly 
perturbed by the small effective coupling. Such weakly perturbed deconfined matter might 
be realized in the core of the neutron stars, or in possible quark stars. However, if the 
temperature of such deconfined system is low enough, the residual attractive force between 
the quarks in the color anti-triplet channel, which is supplied by either one gluon exchange 
(OGE) or the instanton induced interaction, leads to the instability of the Fermi surface. 
This instability causes the non-perturbative generation of the dynamical Majorana mass gap 
near the Fermi surface, which is called "color superconductivity" . Along this line, Bailin and 
Love made an extensive analysis on possible color superconductivities in the quark matter. 4 ) . 

In recent years, there has been great developments in our understanding of the pairing 
phenomenology at high density. The discovery of 11 Color-Flavor Locking" (CFL) state 5 ) at 
extremely high density is one of the examples. CFL state have two types of gap in the chiral 
limit, namely, Ai for the color-flavor singlet excitation, and As for the octet excitations. On 
the other hand, the 2-flavor superconducting (2SC) state have only one gap A, and 4 quarks 
out of 9 obtain this gap. Rich physics contained in the CFL phase, as well as in the 2-flavor 
color superconductivity (2SC) phase, has been revealed 7 \ 

In our present understanding of the high density QCD, pairing patterns other than the 
CFL would hardly take place at high density limit 8 - ) ' 9 **. It is thought that only in the presence 
of the finite strange quark mass m s ^> m Uj( i, the 2SC pairing appears in the low density region 
in the phase diagram. The color-flavor unlocking phase transition from the CFL state to the 
2SC state at zero temperature have been studied using the Nambu-Jona-Lasinio (NJL) type 
model 10 )' n ). This transition occurs when the strange quark mass m s at some fixed chemical 
potential \i gets larger than the critical mass m c s (fi), or equivalently when the chemical 
potential at some fixed m s gets smaller than at some critical chemical potential ii c (m s ). 
Analyses in Refs 10 ^ indicate that a simple kinematical unlocking criterion m^ 2 /4// c ~ 
\A 8 (lic)\m s =o is quite a good guide for the unlocking mass and chemical potential. This 
criterion is based on the observation that the unlocking transition would occur at the point 
where the Fermi momentum for strange quark differs from those for light flavors by the 
smallest gap \A 8 \ for the unperturbed (m s = 0) CFL state. In Ref. 12 \ the unlocking 
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(locking) transition at finite temperature is analysed in the NJL model. 

In the chiral limit, weak coupling analyses of the gap equations reveal the CFL dominance 
over the 2SC at zero temperature and the coincidence of critical temperatures to the QGP 
phase from the 2SC state and from the CFL state. The former is attributed to the subtle 
competition*) between the number of degrees of freedom participating in pairing correlation, 
and the magnitude of the gaps 13 - ) . The latter is due to two facts : One is the rapid vanishing 
of the pairing in the color sextet channel towards the critical temperature T c . The other is 
the disappearance of the nonlinearity of the gap equation in the diquark condensate near 
T c , and which is guaranteed if the transition is of second order. However, in the low, and 
physically interesting density regime, various strong coupling effects such as the participation 
of off-Fermi surface degrees of freedom in the formation of the gap due to the large effective 
coupling constant 14 ^ . Also the pairing in the color symmetric channel may get relevant 
towards low densities, although it is suppressed by one coupling constant relative to the 
anti-trriplet one in the weak coupling regime 13 -* . These strong coupling effects might cause a 
significant modification of these pictures of the CFL dominance at zero temperature and of 
the phase transitions. In addition, critical temperatures to the quark-gluon plasma (QGP) 
phase from the CFL and from the 2SC start to deviate towards low densities, although these 
quantities are shown to be the same in the weak coupling limit 15 **. If this is the case, the 
2SC phase may show up at high temperature region. **) 

For the reasons mentioned above, even in the chiral limit, it would be interesting to study 
the pairing phenomena in dense QCD using some model having strong coupling effects at low 
density side. For this purpose, we examine in this paper, the Schwinger-Dyson (SD) approach 
in the improved ladder approximation to the color superconducting phenomena. 14 - ) ' 18 )~ 20 ) 
We also use the associated Cornwall- Jackiw-Tomboulis (CJT) potential to determine which 
phase of (QGP, 2SC, CFL) is preferred for given density and temperature (/i, T) 7 \ This 
approach, namely the SD equation in the Landau gauge with the improved running coupling 
constant can make the asymptotic behaviour of the mass function at high energy consistent 
with that from the operator product expansion and renormalization group argument in 
QCD 21 \ At the same time, the meson properties on the chiral-broken QCD vacuum, can 
be reasonably understood by the analysis of pseudo scaler Bethe-Salpeter amplitudes 21 )> 22 ). 

*) The condensation energy in the CFL phase is about twice of that in the 2SC phase. If the ratios of 
the 2SC gap and the CFL gaps, A/A\$ get larger than those weak coupling values by a factor \/2 = 1.4 due 
to strong coupling effects, then these two phases may get nearly degenerate with each other in the energy 
density. 

**) In He 3 liquid system, the Anderson-Morel (ABM) state 16 < l shows up in the high temperature and high 
pressure region due to strong coupling effects, even though the Balian-Wertharmer (BW) state 17 ) has larger 
condensation energy in the weak coupling analysis 16 ) . 
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Also, this approach leads to the correct dependence of the gap on the coupling constant 
in the weak coupling perturbative regime. Thus, we believe the improved SD approach is 
reliable for a wide region of the quark chemical potential. Various strong coupling effects, 
say, the large coupling constant, pairing correlation far away from the Fermi surface, and 
that for anti-quark channel, all become relevant for realistic densities in this improved SD 
approach 14 )' 19 ). Also, we investigate the effect of the finite strange quark mass on the pairing 
phenomena, using a simple criterion for the color-flavor unlocking phase transition. 

This paper is organized as the following. Sec. 2 is devoted to the introduction of our 
model and the technical details to derive the gap equations, the CJT effective potentials and 
other physical quantities characterizing the CFL and 2SC states. In Sec. 3, we present our 
numerical results. We discuss the superconducting states at finite T and study the phase 
transitions to normal quark matter in details at /i = 1000 MeV. Also, how the 2SC and 
CFL states at T = and transition temperatures T c are modified towards the lower density 
region by strong coupling effects. At the last of the section, we discuss the phase diagram 
from our model and its dependence on the strange quark mass m s . In Sec. 4, the summary 
of this work and perspective for the future work will be given. 

§2. Gap equations and CJT potentials for 2SC and CFL states 

In this section, we present a general framework to derive the gap equation for the given 
order parameter, focusing on the 2SC and CFL phases. Also, we give various physical 
quantities characterizing the ground state. 

2.1. Gap equations for the 2SC and CFL states 

Nambu-Gor'kov propagators : Here, we introduce two component quark fields \P = (q, g c )* 
in order to express the superconducting gap as a part of the self-energy. In this Nambu- 
Gor'kov base, quark propagators are introduced as 



Here, the subscript c represents the charge conjugation q c = Cq*, and (TO) means the 
ground state expectation value of the time ordered product of the composite operator O. 
We introduce the full self-energy matrix for the quark field in this base as 




^(po,p) = 



E 21 (p ,p) E 22 (p ,p) 



) 
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This self-energy gets dynamically generated by the interaction among quarks and defines 
the asymptotic quark fields. If the interaction is invariant under charge conjugation, the 
relation E 22 (po,p) — —CEn(—po, —p) l C should hold. If the interaction is hermitian so that 
the time evolution of the system is unitary, the other relation E 21 (p ,p) = 70^12 (po, pYlo 
also holds. As a result, four components in S(p) are not independent. 

iS 22 (t,z) = -CiiSni-^-zyyC, S 22 ( Po ,p) = -CS 11 (- Po ,-p) t C, (2-1) 
iS 21 (t,z) = -7o(i£i 2 (-t, -2)) f 7o, S 21 (p ,p) = 7oS , i2(j»o,p) t 7o- (2-2) 

Our main purpose is to investigate the density region where the chiral symmetry is already 
recovered. Therefore, we assume that the diagonal elements of the self-energy, namely, 
Dirac mass function, the deviation of the quark wave function renormalization from 1, and 
the radiative correction to the chemical potential are all zero. We postulate the following 
forms for the off-diagonal elements 4 **' 23 )~ 26 ) ; 

Si 2 (Po,p) = l 5 (A + (p)A_(p ,p) + A4p)A + (p ,p)), (2-3) 
£ 2 i(Po,p) = -l 5 (A-(M-(Po,p) +A+(M+(Po,p))- (2-4) 

A±(p) in Eqs. (2-3) and (2-4) are the projection matrices for positive and negative energy 
states in the massless limit : A±(p) — (1 ± a • p)/2, where a = 707 andj3 = p/\p\. A±(p ,p) 
are the gap functions and are still matrices in the color and flavor spaces. Eigenvalues of A± 
turn out to be gap functions in quasi-quarks (+) and quasi-antiquark (— ) excitations. 75 in 
Eqs. (2-3) and (2-4) guarantees that the pairing occurs in positive parity channel. Although 
the OGE does not distinguish the positive and negative parity channels, the instanton in- 
duced interaction which may be relevant at low densities, prefers the pairing in the positive 
parity channel 27 - ) ' 28 \ The Dirac structures of Eqs. (2-3) and (2-4) are of the most general 
pairing form for the J = + with aligned chirality 23 " ) .*- ) For the moment, we proceed without 
specific ansatz for the color-flavor structure. The inverse of the quark propagator can be 
expressed as 

S (po,p) = „ , v / ; I • (2-5) 
\-E 21 (p ,p) p-iji J 

Here // = /X70 and p is the chemical potential for net quark. By taking the inverse of this 
expression, we get 

c „n , i+ (Po + E +) i>~ (go ~ E -) (oa \ 

£>h(po,p; = +0 — 9 — ^ — jo ~ o — 9 — p« — J7' ^ 2 ' b > 

2p p%-E\- A\ 2p pZ — E 2 — A 2 _ 



*) There is no Cooper instability in the mixed chirality q) channel 29 \ thus the condensate for 

this channel gets suppressed by the factor m q j [i. 
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Su(Po, P ) = -75 U + (p) 2 + A_(p) A- ) , (2-7) 

where we have defined E± = p^ffi and p — \p\, the bare dispersions for quark and anti-qaurk 
in the system with the chemical potential fi. p± are defined as the forward and off-forward 
light cone vectors, p± = (1, ±p), so that jf>± = 7 =F 7 • P- 

Gap equation as a self-consistentcy condition for the self-energy : The gap equations, 
the self-consistency conditions for A±(po,p) up to one loop order, can be expressed as the 
off-diagonal sector of the following Schwinger- Dyson equation in the Nambu-Gor'kov base. 

Efa, P) = W 2 /tJ^ S(q , q)r»(q, p) D%{q Q -p Q ,q-p). (2-8) 
For simplicity, we use the bare vertex instead of the full vertex throughout this paper. 

t a is defined as A /2, where {A a } are Gell-Mann matrices. The vertex r A = — 1^7^ originates 
from the relation qr A q = —q c Cr l A Cq c . Extra minus sign purely comes from the Fermi 
statistics. D^ b u (q — p ,q — p) = 5 ab V ^(qo — p , q — p) is the gluon propagator in the 
medium, and the detailed structure of the function 2^ (go ~ Po,q — p) which we use in this 
paper will be specified latter. 

Through the positive and negative energy projection procedures, we obtain the following 
gap equations, 

d 4 q ( , .t A A + (q ,q)t A 



f d q I 

A±(po,P) = ig 2 Jj^u{c f , u ( y p T ,q + ) 



2 



+C^(P T , g-) ^#% V^(go - Po, q ~ P\ (2-10) 
q - h_ - Zi_ / 

where we have defined c^ v (p,q) by c liu (p±,q±) = -(l/2)tr (A ± (p)^A ± (q)-f 1/ ) 

Color-flavor structure of the gap matrix in the 2-flavor case : In the 2 flavor case, the 
gap matrix in the color-flavor space is specified as 4 ); 

AT c ( Po ,p) = (r 2 x \ 2 )%A ± (p ,p). (2-11) 

Because r 2 A 2 is rank 4 matrix in the color-flavor space, 5 quarks (gu, gd, gs,rs,bs) corre- 
sponding to five eigenvectors with zero eigenvalue of the gap matrix remain gapless. These 
quarks contain the quantum number strangeness or green color, which does not take part in 
pairing. Other 4 modes belonging to the unbroken SU(2) C x SU(2)f doublets [2( r ,b), 2(u,d)] 
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obtain the finite gaps A±. The spectral density matrix p for the 2SC state can be extracted 
from the retarded propagator S§$ c . Its off-diagonal element p\ 2 can be written as 



Pn C (Po,p) = -(r 2 \ 2 )%[y 5 A + (p)A + 27r6(po)5(pl - E 2 + - A\) 

- l5 A_(p)A_27re(p )5( P 2 -E 2 _- Al)] . (2-12) 

The conditions pi — E\ — A±(p , p) 2 = define the quasi-quark and quasi-antiquark disper- 
sions. 

Color-flavor structure of the gap matrix in the 3-flavor case : The color-flavor locked 
state is defined as 5 \ 

if L (p ,p) = (Pi^Ai(p ,p) + (Ps)%A 8 ± (p ,p), (2-13) 
= ^^Ai + (5«5\ - ^5]) Al (2-14) 

The matrix A is full-ranked in the color-flavor space, which results in generating gaps in 
all quasi-quark excitations. *•* In the chiral limit, the CFL state has two-types of gap, corre- 
sponding to two eigenvalues in this matrix, one being A\ and the other being A$. Pus) is 
the projection to the asymptotic quark which transforms as singlet (octet) under the simul- 
taneous inverse rotation SU(3)c+v £ SU(3) C x SU(3)f in the color-flavor mixed space. **) 
The CFL order parameter (2-13) can be divided into states in the color multiplets (3 C , 6 C ) 
as 

A A = -A 1 - -A 8 , A s = -A 1 + -A 8 , (2-16) 
6 3 6 3 v ; 

where N c = 3 is taken. A\ represents self-energy in the color anti-symmetric and flavor 
antisymmetric channel, where the attraction works. However, it turns out that the nonzero 
As gets induced in a non-perturbative way as a result of the self-consistent treatment of the 
gap equation 5 ). One can easily extract the spectral density for off-diagonal propagator, 

PiI L (po,p) = -(P 1 )^ 5 A + A 1+ 2rre(p )6(p 2 - E 2 + - A 2 1+ ) 

+KA_A 1 _2ir6(po)5(pl -E 2 _- A\_) 



*) The CFL states is similar to the Balian-Wertharmer (BW) state 17 ) in the following sense. The BW 
state is characterized by the unit gap matrix in the angular momentum and spin spaces, and has a gap on 
the whole Fermi surface. On the other hand, the CFL state is also a diagonal unit matrix in the color-flavor 
3 C x 3/ space, and this leads non-zero gap for all nine quarks in color and flavor space. The 2SC state is 
much like to Anderson-Morel (ABM) state 16 \ or to polar state, which is not diagonal and isotropic state, 
and has gapless point or line on the Fermi surface. 

**) P$ is a square root of the octet projection matrix (-Ps)ff = <5 ah % — 5f5j/N c . 
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(P 8 )# 75 A h A»27re(p b )<S(pg - ^ - ^ 



8^ 



+ 75 ^_^ 8 _27re(p )5^ ~ E 2 _ - A 2 8 _)\. (2-17) 

The condition p\ — E 2 — Ai^)(p ,p) = defines the dispersions for the SU(3) c +v singlet 
(octet) quasi-quark excitations in the CFL state. 

The Feynman propagator for T ^ can be obtained from above spectral densities, 

7 5 /l + Z\ + (jo ,p) 7 5 yl„zL(j» ,p) 



Sfu(po,p) = 



Po - (\p\ - A 4 ) 2 - A\ + ii] p\- (\p\ + p) 2 - A 2 _+ii] 
+in F (\po\)e(po)pi 2 (po,p; fi). (2-18) 



A and p 12 is color-flavor matrix each for the 2SC and the CFL. n F is the thermal Fermi 
distribution function defined by 

n F (\Po\) = TT ^ ? . (2-19) 

This produces the Pauli-blo eking term in the final gap equation by thermal on-shell quarks 
at finite T, and suppresses the phase space for pairing correlation and brings about the phase 
transition to normal phase at some critical temperature. 

Gap equations for the 2SC and CFL states : Now we simplify the gap equation through 
the vertex decomposition. In the case of the 2SC, we use the following equation, 

8 N + 1 

£ t\Mt A = -^r^A 2 . (2-20) 

A=l ZiVc 

In the case of the CFL, we use the following color-flavor algebra, 

8 1 _ 8 N 2 -1 1 - 

^(Pi)tA = w m J2 i t t A (Ps)t A = ^-(Pi)- Wc (Ps). (2-21) 
These identities can be derived by the Fierz identity : 

8 

£ (AaWAa)^ = -(2/N c )5 ab 5 cd + 25 ad 5 bc . 

A=l 

Using these identities, we obtain 

= K + [A,, A_], jj^A- = K~ [A + , A_] , (2-22) 

for the 2SC, and 

2N AN 

' A l+ = -K + [A 8+ ,A 8 -], 2N C A 8+ + -^-A 1+ = -K + [A 1+ ,A 1 _] , (2-23) 



N 2_l ^ L -o-r,-o-j, -o-ot ■ N 2_ 1 

2N 4N 

= ~ K ~ A ^ 2N ^~ + W=T Al ~ = ~ K ~ [Al+ > Al - ] ■ (2 ' 24) 
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for the CFL. We define the kernel integral for OGE as: 

K a [A + ,A_] = -ig 2 j ^{s$(q;A + )c ia ,(p*,q + ) 



+S F (q;A_)c^(pa,q-)^V^(q -p ,q-p), (2-25) 



with definition 



S$(q;A ± ) 



Qo-(\q\T P) 2 - A± + iv 

+2mn F (\q \)5 \q 2 - (\q\ =F /i) 2 - A%] . (2-26) 



Index a takes + and — , and we define a = —a. These gap equations are the generalization 
of those obtained in Ref. 5 ** to the non-local interaction in arbitrary N c with the antiquark 
gap. 

Gluon propagator in the quark Fermi liquid : We use the quasi-static approximation of 
the hard dense loop propagator in the Landau gauge, as it is taken in Refs. 14 )> 24 )> 30 ). 

2V(Po,p) = -A T (Po,P)P^ - MPo,p)P^, (2-27) 
with the transverse and longitudinal projection matrices, 

Pij = °ij ~ PiPji P(H = PiO = ^00 = 0? Pfiv = ~9{i.v + ~n — P^ u , (2-28) 

Po ~ P 

and the corresponding amplitudes 

A T {p ,p) = 2 . . Q , T7 2 TTTi | / , A L (p ,p) = 1 (2-29) 

p 2 + t9(p-p )(rnn 2 D /A)\p \/p p 2 + m 2 D 

m 2 D = Nfg 2 /i 2 /2ir 2 is the Debye screening mass in the quark-gluon plasma. In magnetic 
sector, we have to take into account the subleading correction from the dynamical sector 
up to linear in p /p in order to avoid the collinear singularity. *) We have neglected the 
effect of the temperature on the polarization because our main interest lies in the region 

Z\(p)| T =o< pi- 
Gap equation for quasi-particle on the mass shell : Now we simplify the gap equation 
by integrating out the frequency as in Ref. 14 ^ because our main interest lies in the gap 
of quasi-quarks on the mass shell. In the kernel K a of the gap equations, we have three 
contributions: quark pole, plasmon pole and blanch cut in the gluon continuum. It has 

*) The quasi-static approximation can be justified if the ratio of the obtained gap to the Fermi surface 
(A/ /x) is small enough to neglect (Z\//x) 2 . The quasi-static approximation is going well in our model, because 
it will turn out Aj \x ~ 1/10 even at the lowest density. 
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been shown in Ref. 24 ** that only the quark pole has a significant effect on gap equations at 
least in the weak coupling limit A(p)/fi — > at T = 0. Here, we assume this property 
holds even at finite temperatures and/or in the strong coupling regime. The resulting gap 
equations for the 2SC and CFL states under quark pole dominance take the same forms as 
Eqs. (2-22)~(2-24) under replacement of the kernels K a by : 



K± [p; At] 



dq T) f n „\ A+ ( q 



+ 



(2tt) 3 
dq 



(27T) 



p±-(p, q) 



2e + {q 
A_(q 



1 - 2n f 



2e_(q) 



1 -2n F 



(9) 



(2-30) 



where e±(p) = \J{pT /-0 2 + A±(p) 2 , and A±(p) are gaps on the quasi-quark and quasi- 
antiquark on the mass shell. The effective interaction T> a/3 (q,p) is constructed by the fol- 
lowing electric and magnetic interactions, 



V a(S (p,q)=V^(p,q) + 2V™(p,q), 



(2-31) 



with 



+± 



[p,q) 



V%{p,q, 

TRe 
4 



9 2 (l±p-q)/2 
(p-q) 2 + m 2 D ' 

g 2 ((p -q) 2 ±(p-p-q)(q-g- p)) 

(p - q) A + i{nm 2 D /A)\e a (p) - e±{q)\\p - q\ 

g 2 ((p - qf ± (p - p ■ q)(q - q ■ p)) 



iO{g 4 



(p - q) 4 + i(irm 2 D /4)\e a (p) + e±(q)\\p - q\ 

The Fermi distribution function n F in Eq. (2-30) represents Pauli blocking by 
excited quasi-quark and quasi- ant iquarks. After integration over angular variables, 
the following kernels. 



thermally 
we obtain 



K^p- A ± 



+ 



9 



16tt 2 



f°° q 
/ dq— tanh 
Jo k 



£+\ A + (q) \ {p±q) 2 + m 2 D ( (p + q) 2 + m 2 D 
2Tj2e + {q)\ " 2pq ' \{p - q) 2 + m 2 D 



l ln ( {p + qf + M\e ± {p)-t + {q)) 2 (p + g) 6 + M\e ± {p) + e + (g)) 2 
3 n l ( p _ g )6 + M 4 (e±(p) _ e+(g))2 ( p _ g )6 + M 4 (e±(p) + e+(g))2 



9' 



+ — 2 f^tanhf^^44^ 



(p-q) 2 + ml ( (p + q) 



2pq 



(p - q) 2 + m 2 D 



+- 



In I 



(p + g) 6 + M 4 (e±(p) - e_(g)) 2 (p + g) 6 + M 4 (e±(p) + e_(g)) ; 



(p - g) 6 + M 4 (e±(p) - e_(g)) 2 (p - g) 6 + M 4 (e±(p) + e_(g)) ; 
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where we have defined M 2 = nm 2 D /4. We have neglected the parts with higher angular 
momenta in the gap function A(q).*^ Following Ref. 14 \ we replace the coupling g 2 with the 
momentum dependent effective one g 2 (p, q) by the Higashijima-Miransky prescription 21 \ In 
the improved ladder approximation, g(p, q) 2 is taken to be 

16tt 2 1 



9 ip, q) = 



Pmax = max(p,g), 



(2-32) 



(3 ln((pLx+^)M 2 ) : 

where @ = (11 N c — 2iVy)/3, pi plays a role of a phenomenological infra-red regulator. We 
adopt A = 400 MeV and pi = 1.5 A 221 ) for numerical calculations. 



2.2. CJT effective potentials for finite temperature 

The CJT effective potential 31 ) (The Ward-Lattinger effective potential in the case of 
non-relativistic fermion) for fermion propagator is given by, 



r [S, p] = \ {TrLog [S] - Tr [SS 1 ] + V[S]}. 



(2-33) 



Here, "TrLog" should be interpreted in the functional sense. In the momentum space, 
S^(po,p) = S 1 (p ,p) - H(po,p) with Sq\po,p) = diag.(#f + ^,#f -^/). The overall factor 
1/2 in the Eq. (2-33) should be introduced to eliminate the artificial degrees of freedom 
introduced by the Nambu-Gor'kov 2-component fermion. \i represents explicit /i dependence 
through the bare propagator and the gluon polarization. V[S] is a potential functional which 
contains Feynman diagrams of the skeleton self-energy and one quark propagator S. This 
is a generating functional of the proper self-energy, namely, the functional derivative of the 
potential is the self-energy, 5V[S]/5S = S[S}. Now we take one loop approximation for 
E[S], symbolically, written as E {1) [S] = g 2 r SrD.**> In this case, r[S] reduces to 



r (2) ~s. fl ]- ] -\ II,,, S 



Tr 



SS, 1 



+ ^Tr \SS^ [S] 
2 



(2-34) 



The saddle point approximation leads to the one loop gap equation : ST^/SS = S 1 — 



S 1 + 2^ [S] = 0. We write the solution of this equation S(g; /j) and S 1 — S 
At the stationary point, we can simplify the effective potential: 



r (2) 



S(g;p),p\ = - \ TrLog 



-Tr 



SU 



E{g;p). 



(2-35) 



This is a useful formula in which the coupling and gluon propagator do not appear explic- 
itly, but are hidden in S and E. Using CJT method described above, the thermodynamic 

*) The gap function A(q) can be expanded as A(q) = A(q) A 1 (q) Pi (cos 8). We assume A 1 < A, 

because the isotropic gap function can activate the whole Fermi surface where the infrared Cooper singularity 
lies. 

*) r and r means the bare and full vertices in the Nambu-Gor'kov base. 
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TrLog[S]-Tr[sV] V (2) [5] 

Fig. 1. Schematic Feynman graph of 2-loop CJT potential functional. S-derivative on this poten- 
tial is removing one solid line, which yields one loop gap equation. 



potential (condensation energy) density relative to the normal Fermi gas up to 2-loop order 
is given by 



6f2(n,T) 



■TrLog 



-i T r 
4 



[2-36) 



PV 12 

Note that this expression is only valid at the stationary point which can be determined 
by the gap equation for S. After an involved but straightforward algebra and Matsubara 
summation, we reach the following expression for the 2SC phase, 



sa 



2SC 



2 N C -1 

3 N r . 



-T 



dq 



In 



1 

+ 2 



dq (A 2 



(2tt) 3 I 2e. 



tanh — 



£+ 
2T 



cosh(e+/2T) 
cosh(|£+|/2T) 

A 2 



+ In 



cosh(e_/2T) 
cosh(£_/2T) 



2e_ 



tanh ( — 
2T 



(2-37) 



Roughly speaking, the first term represents the thermodynamic energy gain by the dynamical 
quark loop in the presence of the condensate, while the second term represents the energy 
cost to put the condensate in the Fermi gas. 2(N C — 1)/3N C in front of right hand side of the 
above equation indicates that the flavors joining the pairing are 2 out of 3, while the active 
colors are N c — 1 out of iV c . 

In the CFL case, the algebra is more complicated but a straightforward calculation yields 
the following result, 



5a 



CFL 



6N n 



-T 



-T 



1 

+ 2 



+ - 



dq 

(2^F 

dq 
(2^)3' 

dq 

dq 



-Lin 

N 2 



N 2 



In 



cosh(ei+/2T) 
cosh(|£+|/2T) 
cosh(ei_/2T) 
cosh(£_/2T) 



N 2 -l 



N 2 



In 



+ 



N 2 

c 



N 2 



In 



cosh(es+/2T) 
cosh(|£+|/2T) 
cosh(e 8 _/2T) 



N 2 2e 1+ 
1 A\_ 



2 J (2tt) 3 \N 2 2e^ 



tanh 



tanh 



(- 

\2T 

(— 

\2T 



+ 



N? - 1 A 



8+ 



TV 2 

c 

N 2 - 



2e 



8+ 



1^1- 



cosh(.EL/2T) 



tanh ( 

2T 



N 2 2e 8 „ 



tanh 



2T 



(2-38) 



Eqs. (2-38) and (2-37) take simple forms 5f2 CFL (/j,T)/i2 4 = f{A l /n) + {N c N f -l)f{A s /^) = 
f(A 1 / / j,) + 8/(zVm) and <^ 2 sc(m, T)//x 4 = 2(N C - l)f{A/n) = Af(A/fi). Therefore, if the 
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gaps A, A\ and Z\§ take all the same value, then the condensation energy for the CFL phase 
is 9/4 times larger than that for the 2SC phase. The comparison of the (u, d) + s 2SC matter 
with the (u, d, s) CFL phase is not quite obvious till one adopts a model and solves the gap 
equations and the thermodynamic potential. 



2.3. Occupation numbers and Correlation functions 

Here we give expressions for various quantities in the ground state. 

Quark number density : Occupation number characterizes how much the Fermi sphere 
is distorted by the gap. We can extract the occupation numbers from the frequency integral 
of the diagonal propagator Su(qo, q), 



r,b/ \ ' -^m 
n+ (g) = 2-27 



^tanhfeV n r ' b (q) = - 
2e+ \2T ) ' w 2 



1 , /e_\ 

tanh — - 

2e_ \2TJ 



(2-39) 



n 9 + (q) = (l-n F (\E + \)9(-E + )+n F (\E + \)9(E + ), n 9 _ (q) = n F (E). (240) 

n + (q) and n_(g) are distribution functions for quark and antiquark. n F (\E + \)9(—E + ), 
n F (\E + \)9(E + ) and n F (E_) represent the thermal hole, the thermal quark and the anti- 
quark, respectively. Net quark density for the 2SC in some chemical potential ji can be 
obtained by 



M = 4 / 



dq 

(2^ 



(N c - l)(nt(q) - n r *{q)) + (n 9 + (q) - n 9 _(q)) 



(241) 



The chemical potential has to be determined by the number conservation p = p(/x), which 
also can be obtained by the derivative of the pressure — i?(p). 
Similarly, we can find occupation number in the CFL phase as 



n 1± (q) = ^tanh ( , n 8 ±(o) = - 

W 2 2ei± \2T J 8±w; 2 



— — tanh I j . 

2e 8± \2TJ 



(242) 



rii and rig are the color-flavor singlet and octet distribution, which can be defined using the 
annihilation operator for quark, a{p) and that for antiquark b{p) as 

(af(p)a%p)) = (POjfm+Cp) + (P 8 )-N+(P), 



^(P)6J(P)> = (A)?rn_(p) + (P^ffns-Cp). 

Using these octet and singlet densities, we find by putting b = a and j — i in the above 
equation, 



\ ab 



ni±(q) - n 8 ±(q) 



(243) 
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nf(p) is the probability distribution finding a quark (+) or an antiquark (— ) with color a 
and flavor % having the momentum p. We can see the fact that the quark density is almost 
determined by the octet density. The total quark numbers for quark and antiquark read 



dq 



(N C N F - l)(n 8+ (q) - ra 8 _(g)) + (n 1+ (q) - m_(g)) 



(2-44) 



Correlation functions : What is important to characterize a superconducting phase is a 
energy gap and the correlation functions which describes the off-diagonal long range order 
in the system. This quantities are also called the anomalous density because it is related to 
the anomalous propagators (Tqq). In the case of the 2SC, 



0±(g) = tanh ( ^ J 



v2T/ 2e± 

Here, the correlation functions <p±(p) are defined by 



1 - 2n F {e±) 



2e±' 



(2-45) 



«(PK(-P)> = (r 2 A 2 )^ + (9), (bttp)b b j(-p)) = (r 2 A 2 )5^_(9). 

To avoid an involved notation, we have omitted the helicity indices. As is postulated for the 
ansatz of the 2SC condensate, the correlation exists only in the color-flavor anti-symmetric 
channel. The minus sign in front of np in Eq. (2-45) shows the effect that the thermally 
excited quasi-particles reduce the strength of the correlation. 

In the same way as in the 2SC case, we easily derive the correlation functions in the CFL 
phase, The anomalous densities in the singlet and octet channels become 



</>i±(q) = tanh 



(—^ — ^, (j>s±(q) = tanh ^ " 1 



\2T J 2e 1± 



2TJ 2e 8± 



(2-46) 



with definition 



(a?(p)aj(-p)> = (Pi)i<h+{q) + (A)#6h-(?), 

<6?(p)6J(-p)> = (Pi)^;_(?) + (Ps)^;m- 



By using the locking relation ej^e 1 ^ = (N c — l)(Pi)"j > — (Pg)"]', the color 3 C flavor antisymmetric 
correlation (/>a(q) and the color 6 C and flavor symmetric correlation 4>s{q) are extracted as 



We would see 05 disappears rapidly near the critical temperature compared with 0^, due to 
anti-triplet dominance in the critical region. 
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3. Numerical results 



In this section, we present our numerical results. First, we focus on how the temperature 
affects the superconducting gap, and other physical quantities at a high chemical potential 
li = 1000 MeV. We discuss various features of the phase transition. Then, we see how 
these pictures of the phase transition would be modified towards low, physically interesting 
densities. Finally, we draw the QCD phase diagram for the region \x > 400 MeV. We estimate 
the effect of the effective strange quark mass m s on the QCD phase diagram using a simple 
unlocking criterion, and study how the QCD phase diagram changes with a variation of m s . 

3.1. Phase transitions to QGP phase at /i — 1000 MeV 

Gap functions and condensation energies : We show the numerical solutions for the gap 
equations Eqs. (2-22)~(2-24) at a fixed chemical potential /i = 1000 MeV corresponding to 
about 80 times the normal nuclear density p — 0.17 fm~ 3 . Fig. 2 shows the gap function 
A±(p) for the 2SC state (a), and the singlet and octet gap functions (Ai±(p), A 8 ±(p)) for the 
CFL state (b). At first sight, we can see the shapes of these gap functions are all similar as a 
function of momentum except for the magnitude |Z\ 8 (jo)| < < |Al(p)| and the relative 

sign between the singlet and octet gaps. Secondly, we observe the temperature uniformly 
reduces the magnitudes of the gap functions towards the critical temperature which can be 
read ~ 15 MeV. 

In order to see characteristics of phase transitions in more detail, we show in Fig. 3, the 
temperature dependence of the gaps at Fermi surface (a) and that of the CJT condensation 
energies for the 2SC and CFL states (b). From these figures, we can see that the phase 
transitions to the QGP phase both from the 2SC state and from the CFL phase are of 2nd 
order just like in the BCS theory with contact interaction : The derivatives of the gaps seem 
to diverge towards the T c in (a), and thermodynamic bulk quantities in superconducting 
phases are connecting continuously to those of normal QGP phase (b). Another thing to be 
noticed is that the critical temperatures to the QGP phase are the same both for the CFL 
and for the 2SC. We will discuss this property in more detail latter in this section and see 
that this holds at least in the weak coupling limit. 

The condensation energy for the CFL state is always about 1.9~ 2 times larger than 
that for the 2SC in spite of the fact that the octet gaps determining the thermodynamics in 
the CFL takes smaller values than the 2SC gaps. This is consistent with the weak coupling 
analytical estimate 13 ^ : The CFL phase is easier to excite quasi-quarks, but the bulk energy 
is lower in the CFL than in the 2SC. This reversal between the condensation energy and the 
magnitude of gap is clearly a direct consequence of the fact that the number of the active 
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Fig. 2. Gap functions at fi = 1000 MeV and various temperatures. We have drawn the gap for 
antiquarks at negative momentum region, namely, (— |fc|, Z\_(|A;|)) are plotted at left half of the 
figure, (a) Gap function for the 2SC state, (b) Gap functions for the CFL state. 
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Fig. 3. CFL vs. 2SC. 

(a) Gaps at the Fermi surface vs. temperature, 
quark for both phases at fi = 1000 MeV. 



(b) 
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Thermodynamic potential and energy gain per 



freedom taking part in pairing is 9/4 times larger in the CFL phase than in the 2SC phase*) 



*> The condensation energy for the CFL is a function of A\ + 8Ag ~ 12Zi|, and that for the 2SC is a 
function of 4A 2 . Thus, the thermodynamic enhancement factor 9/4 is a little bit underestimated, but is 
enough for explaining the CFL dominance at T = 0. 
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Fig. 4. (a) Correlation functions for the 2SC at finite temperature, (b) The coherence length 
and the Pippard length in the 2SC state. Black line show the thermal wave length of the 
un-perturbed Fermi gas ~ 1/2T. 



It can be said that the CFL state is relatively more fragile.*- 1 than the 2SC state due to 
smaller gap in the color-flavor octet channel. Reflecting this, the critical temperatures to 
the QGP phase from both states are almost the same although the condensation energy for 
the CFL at T = is 1.9 times larger than that for the 2SC state. 

Cooper pair size at fi = 1000 MeV : In Fig. 4(a), the correlation functions for the 2SC 
at various temperature are shown. These quantities characterize the internal structure of 
pairs. Although the magnitude of the correlation gets significantly reduced by temperature, 
the width of the function seems not so much affected. This implies that thermally excited 
quasi-qaurks destroy the phase coherence of the system, but does not change the internal 
structure of pairs. This fact gives us a picture of the phase transition : The phase transition 
occurs by the reduction of the number of coherent pairs, not by the dissociation of them. 
In order to make sure of this point, we have drawn various length and those dependence 
on temperature in Fig. 4(b). As is speculated above, the coherence length £ c denoted by 
black points are not much affected by temperatures. £ c is defined by the root mean square 
radius of the Cooper pair with wave function <p + . Black line shows 1/2T, which is related 
to the thermal de Broglie wavelength of the naive Fermi gas : At = \j \j (E\)t ~ 1/2T. 
Only the interaction with the length scale < Xt can cause quantum corrections to the naive 

*) Because the quasi-quarks contribution to the thermodynamic quantities in the CFL can be expressed 
as the integral of 8e~ e8 ^ T + e _€l / T , the color-flavor singlet quarks with the large gap do not give a signif- 
icant thermal effect for T < T c . On the other hand, in the 2SC state, thermal quark contribution to the 
thermodynamic quantities can be written as the integral of 4e _e / T . Thus, the CFL state is more sensitive 
to temperature than the 2SC state because of octet quasi-quarks with the smallest gap. 
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Fermi gas. Phase transition to the QGP phase takes place at these 2 length scale becomes 
comparable, because the quantum correlation with length £ c gets gradually destroyed by 
thermal disturbance with length At- 

Why critical temperatures coincides? : Now, we return to the problem why the critical 
temperatures to the QGP phase coincides. It is related to the disappearance of nonlinear 
effects near critical temperature T c . To see this, we rewrite the CFL gap Eqs. (2-23)~(2-24) 
in terms of the symmetric and anti-symmetric self-energy (A a, As) in Eq. (2-16). If we 
neglect the Ag term in the loop integral in Eqs. (2-23) and (2-24) (ansatz of the anti-triplet 
dominance \A$\ <C |4i|), and also the contribution from antiquaries which can be justified 
only in the weak coupling region, equations for Aa and As cast into the following form: 

4TV 2 

-A A = (N + 1)K + [A A ]+K + [(N-1)A A ], (3-1) 



TV + 1 

' " ^ A s = (TV - l)K + [A A ] - K + [(TV - I) A A ] . (3-2) 



TV- 1 

If we neglect the nonlinear terms in Aa, we obtain 



47V 2 47V 2 

-A A = 2NK + [A A ] , ——A s = 0. (3-3) 



TV + 1 TV — 1 

These two "gap equations" are consistent with the first ansatz of omitting As in the loop 
integral. Discarding nonlinear terms in A a is justified near the critical temperature of the 
2nd order transition. To make sure of the validity of this argument, we check the ansatz of 
the color anti-triplet dominance near the critical temperature. Fig. 5(a) shows correlation 
functions in the anti-triplet and sextet channels. The sextet correlation disappears rapidly 
towards the critical point compared to the anti-triplet correlation. Similarly, Fig. 5(b) shows 
the temperature dependence of the ratio —A 8 /Ai. This quantity approaches 1/2 near the 
critical temperature, which directly means the anti-triplet dominance As/A A — > with 
T -> T c . The convergence -4/4 -> 1/2 seems linear in t = (T c - T)/T c , therefore, if 
A A (t) = —Ag(t) ~ t 1 / 2 near T c just as the scaling in the mean field approximation, then As 
approaches as A s ~ t z l 2 towards T c . Namely, As disappears rapidly towards T c . 

3.2. Gaps, condensation energies and critical temperatures for ji = 1000 ~ 400 MeV 

Now we turn to the axis of chemical potential in the phase diagram, and investigate how 
the physical quantities studied above change, and how the picture of phase transition gets 
modified as one goes to lower density (/i < 1000 MeV) region. 

Gaps and energy gains per quark at T = : We show the chemical potential dependence 
of the 2SC gap A and the CFL octet gap \A 8 \ in Fig. 6(a). The ordering \A 8 \ < A is true 
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Fig. 5. Disappearance of nonlinearity near critical temperature, (a) The temperature dependence 
of correlation functions in the anti-triplet and sextet channels for the CFL state, (b) The 
temperature dependence of the ratio —A%/A\. 



for all region 300 MeV< [l < 1000 MeV. Our gaps have several features which are quite 
different from those obtained in Ref. 5 ). 

(1) The gaps in Ref. 5 ) disappear for /x > 800 MeV. Because the NJL model used in the 
reference is an effective theory originally introduced to describe the chiral symmetry break- 
ing in the QCD vacuum, it has a cut-off scale A ~ 1 GeV. Our gaps are also decreasing 
functions of /x, but do not go to zero even at high density. Furthermore, our gaps will turn 
into increasing functions of /x at some higher chemical potential mu > 1000 MeV, because 
our model coincide with other perturbative SD approaches 32 )> 25 )> 26 )i 24 ) ) which all predict 
increasing gaps in the weak coupling regime. 

(2) The gaps obtained in Ref. 5 ) go to zero towards /x — > 0, while our gaps do not. Gap 
equation at /x = is the same as the chiral gap equation under the replacement the NJL 
4-Fermi coupling K — > K/2, if we neglect the small admixture of the pairing correlation 
in the color and flavor sextet channels. The coupling K has been tuned to reproduce the 
chiral gap 400 MeV at zero density. Therefore, the disappearance of the superconducting 
gap at /x = means that K/2 < K c < K, where K c is the critical coupling for the dynamical 
generation of the chiral mass gap. On the other hand, our gaps increase as /x — > and seem to 
converge to non-zero finite values. It can be said that our model describes a strong coupling 
superconductor at low densities in the sense that the large effective coupling constant leads 
non-zero gaps even in the absence of the Fermi surface 14 ^. Of course, if we had included the 
diagonal self-enrgy part which is relevant to the chiral symmetry breaking, then one would 
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Fig. 6. 2SC vs. CFL at various densities, (a) Chemical potential dependence of the octet gap and 
the 2SC gap. (b) Density dependence of condensation energy gain per quark both for the CFL 
and 2SC states. 



encounter the phase transition to chiral broken phase at some critical chemical potential /j, c , 
and the superconducting gap would disappear at zero density due to the reduction of the 
state density at Fermi surface owing to the generation of mass gap 33 - 1-35 ''.*-' 

In Fig. 6(b), the chemical potential dependence of the energy gain per quark relative to the 
naive Fermi gas is shown. First, the CFL state persists against the 2SC for all density region, 
and its energy gain is about 1.9 times larger than that in the 2SC. Second, this quantity 
also shows an increasing behaviour as one goes to lower density region. Nevertheless, it is 
small in order of ~ lMeV even at fi = 300MeV. This fact is attributed to the fact that 
particle energy gain is proportional to ~ A(A/fi), namely, the pairing effect on the energy 
gain per quark is suppressed by factor Aj '//. This means that only Aj jj, of total quarks near 
the Fermi surface get the energy gain —A. However, Z\//i~ 1/10 at fi = 300 MeV is quite 
larger than the case of the BCS weak coupling superconductor, where typically this ratio 
takes the value ~ 1/1000. It can be said that the color superconductor which might be 
realized in the neutron star core is high T c one of strong coupling. 

Critical temperatures : Now we discuss the chemical potential dependence of the critical 
temperature. In Fig. 7 (a), critical temperatures for the 2SC and CFL states as a function 
of the quark chemical potential are shown. The dependence of T c on the chemical potential 
is similar to those of gaps in Fig. 6(a). The strong coupling effects also exist and lift up 
critical temperatures at lower densities. Even then, the critical temperature at the lowest 
density is ~ 20 MeV. Although the behaviour of the critical temperature as a function of 

*' However, there is a model in which the diquark gap persists against the sudden jump of the chiral 
gap 36 ) . 
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Fig. 7. (a) Critical temperatures extracted by fitting the data in the from A(T) = a(T c — TY 

chemical potential (b) The deviations of the ratio T c (u)/A(fi,T = 0) from its weak coupling 
universal value £, vs. A/fJ,, a variable for the weak coupling expansion. 



jj, is qualitatively similar to that obtained from the improved SD approach for m s = oo 
in the Ref. 2 °), the magnitude of the critical temperature is much smaller in our case than 
~ 150 MeV in the reference. The correct treatment of the Landau damping in the magnetic 
gluon would make this difference smaller. Fig. 7(b) shows the ratio T c (fi)/A(fjL,T = 0) as a 
function of A(fi,T = 0)//i. This ratio takes the coupling independent universal value in the 
weak coupling limit Aj\i — > in the BCS model with a contact interaction. We can see that 
this ratio approaches to the weak coupling value towards A/fi —>■ 0. However our interested 
region is far away from the region where the convergence to the weak coupling values can be 
seen. Again, we notice the importance of the strong coupling effects. 

Even in the presence of these strong coupling effects at low densities, the ratio of the 
critical temperatures for the CFL and 2SC states do not deviate from its weak coupling 
value. Namely, the critical temperatures for the 2SC state coincides with that for the CFL 
state for all region of chemical potential within numerical errors. 

Simple estimate of the unlocking transition for T = : Here we estimate the color-flavor 
unlocking transition at T = using a simple kinematical criterion, which turned out to 
be quite a good guide for unlocking transition in the NJL-like model analyses 10 - ) ' In 
Fig. 8, we have shown the critical strange quark mass in a kinematical criterion which can 
be extracted from our result for the unperturbed octet gap as m^(fi) = 2^J /iZ\ 8 (/i)| ms=0 ,. 
Suppose that the effective strange quark mass m s (fi) is decreasing function of the chemical 
potential \i. If m s (300 MeV) > m c s (n) and m s (1000 MeV) < m^(/i), then m s (//) intersects 
the curve of m c s {n) at some point /i c locating between 300 MeV < \i < 1000 MeV. In the 
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Fig. 8. The chemical potential dependence of the critical strange quark mass from a simple crite- 
rion for the unlocking transition with Agfa) obtained by solving the CFL gap equations. 

simplified (kinematical) picture, the unlocking transition occurs when the chemical potential 
\x gets smaller than this /i c , where the deviation of the Fermi momentum for strange quark 
from those for light flavors, p s F — pp d = \j\i 2 — m s (fi) 2 — fi ~ m s (fi) 2 /2/x, exceeds twice of 
the smallest gap \A 8 \ in the CFL state. Here, we had better mention that if we take into 
account the dynamical effect of the strange quark mass on the formation of gaps, then the 
gap of quarks belonging to color-flavor octet splits into 4 different values, and the smallest 
one Al at the critical point sets the true unlocking critical mass m^(true) = 2-Jfj.A^. In the 
Ref. u \ the deviation of the smallest gap from the unperturbed octet gap As is about 20% 
at fi — 400 MeV. Thus using the simple kinematical criterion, we might overestimate the 
critical mass about 100(1 - VoI)/VoI ~ 10% of the real value for /i = 400 MeV. 

In the following, we assume that this criterion is useful even for finite temperature, and 
we put the further ansatz that the effective mass of the strange quark at chemical potential 
fi can be approximated to constant m s (fi) = m s for region 300 MeV < /i < 1000 MeV. 
Instead of this simplified ansatz, we treat m s as a parameter, change it by hand, and discuss 
the qualitative dependence of the phase diagram on m s . 

3.3. Phase diagram and its m s - dependence 

Here, we draw a phase diagram expected to be realized from fi ~ 300 MeV to fi = 
1000 MeV for various values of the strange quark mass. As already noted above, we consider 
constant effective strange masses without /i-dependence. Fig. 9 shows our result of the QCD 
phase diagrams for four values of the strange quark mass m s = 0, 150, 200 and 250 MeV. 

1. The chiral limit m s — : In the chiral limit m s = 0, there is only one line (black bold dot) 
which divides the (/i, T) plane into the QGP and CFL phases by the 2nd order transition. 
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This degenerates in double lines, namely, one separating the QGP and the CFL, and the 
other dividing the QGP and the 2SC. However as we discussed above, the 2SC phase does 
not appears for m s = 0, because the condensation energy of the CFL is about twice larger 
than that of the 2SC for all region T < T c (fi). 

2. m s = 150 MeV : We study what happens if we switch on the strange quark mass 
m s = 150 MeV. We estimate this using the kinematical criterion descried at the end of 
the previous section. New line denoted by triangle-line appears in the phase diagram. We 
call this color-flavor unlocking critical line. This unlocking line divides the superconducting 
region T < T c (/x) into the CFL and 2SC phases by 1st order transition. Above this unlocking 
critical line, the system is unlocked into the 2SC state, and below it, the system stays in 
the CFL phase. On the unlocking critical line, the pressure for the CFL coincides with that 
for the 2SC. The strength of the 1st order phase transition gets weaker and weaker towards 
/i — > oo, because the gap between energy density for the 2SC phase and that for the CFL 
phase gets smaller as one going near the unlocking line. The critical end point at which the 
1st order transition terminates is located at fi = oo as long as we use a kinematical criterion, 
but the dynamical effect of the strange quark mass on the pairing may bring this point to the 
finite density point. Indeed, in the NJL model analysis 12 **, the tricritical point appears at 
fi ~ 520 MeV, although the physical reason of the appearance of 2nd order phase transition 
has not yet been understood. Furthermore, because this tricritical point is quite close to the 
cut-off A ~ 600 MeV, it might be an artifact of such an effective model with a cut-off scale. 

3. m s = 200 MeV : Now, we change the strange quark mass m s = 150 MeV to m s = 
200 MeV. Then the unlocking critical line shifts from the "triangle-line" to the "square-line" . 
The room for the 2SC realization is enlarged, pushing the CFL area to higher density and 
lower temperature region. Here we assume the chiral transition occurs near \i = 400 MeV, 
then at zero temperature, the system is always chirally broken by the qq condensate for 
ii < 400 MeV, and by the CFL type qq condensate for fx > 400 MeV. In this situation 
m s < 200 MeV, quark-hadron continuity might possibly realize. 

4- m s = 250 MeV : Finally, the "cross-line" represents the unlocking critical line for m s = 
250 MeV. In this case, the 2SC might intrude into the opening region between the chirally 
broken vacuum phase and the CFL phase in T = section in the QCD phase diagram. 
Quark-hadron continuity does not occur for m s > 250 MeV in our model. The critical mass 
for the quark-hadron continuity lies in the interval 200 < mf nt < 250 MeV in our model. 
The dynamical effect of the strange quark mass on the pairing gap would bring this mass 
to lower densities. The true critical mass m^ ont from our model would be a little bit smaller 
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Fig. 9. The QCD phase diagram for various values of the strange quark mass. The bold-dot- 
line represents the phase boundary between 2SC and QGP separated by the 2nd order phase 
transition. Triangle-line, square-line and cross-line are showing the gradual change of the 
unlocking critical line when we shift the strange quark mass as m s = 150, 200 and 250 MeV 
by hand. 



We have applied the Schwinger-Dyson equation in the improved ladder approximation 
to the color superconducting phenomena in high density quark matter. Special attention 
has been given to the competition between the 2SC phase and the CFL phase at finite 
temperature, and to the effect of the strange quark mass on the QCD phase diagram. 

Sec. 2 was devoted to the technical details to obtain the gap equations and the CJT 
potentials for CFL and for 2SC from the Nambu-Gor'kov Schwinger-Dyson approach. In 
Sec. 3, we studied which ordering of (2SC, CFL) is favored at relatively low density /J = 
300 ~ 1000 MeV at finite temperature. We have drawn the picture of the phase transition 
by investigating the \ijT dependence of the gap function, the correlation function and the 
coherent length, as well as bulk quantities such as pressure. We summarize our main results 
in this paper below. 

(1) The phase transitions to the normal phase both from the 2SC and from the CFL have 
shown almost the same characteristic behaviour as in the BCS weak coupling superconductor 



than ~ 250 MeV obtained from the N JL model in Ref . n > . 



4. Summary and Discussion 



24 



with contact interaction. The critical temperatures for the CFL and the 2SC coincides with 
each other due to the color anti-triplet dominance, despite the fact that the condensation 
energy density is almost twice larger in the CFL than in the 2SC. This facts are consistent 
with weak coupling estimates, and are indicating that the CFL state is more fragile than 
the 2SC state against thermal fluctuations, due to the smaller energy gap in the color-flavor 
octet channel. Phase transition occurs as a consequence of the reduction of the coherent 
Cooper pairs, but not due to the dissociation of those. 

(2) The strong coupling effects get larger towards the lower density side. Especially, this 
was observed in the violation of the universal relation between the gap at T = and the 
critical temperature. This value deviates from the weak coupling analytical value ~ 15% at 
(j, = 300 MeV. This deviation comes from strong coupling effects, the large off-Fermi surface 
effect including the pairing in the antiquark channel owing to a larger coupling at lower 
densities. We find that physically interesting densities lie far away from the region where 
this ratio seems to converge to the universal value. Despite these strong coupling effects, 
the ratio of the physical quantities in the CFL state and in the 2SC state, does not so much 
differ from the weak coupling value even in the strong coupling region. 

(3) The QCD phase diagram and the effect of the strange quark mass m s on it are stud- 
ied within a simple unlocking criterion. In the chiral limit, there is no room for the 2SC 
realization, but the finite strange quark mass makes a unlocking critical line separating the 
superconducting (fx, T) area into the 2SC phase and the CFL phase, and this line push the 
CFL region to higher density and lower temperature region. The critical end point on which 
the 1st order unlocking line terminates is located at fi = oo as long as a simple kinematical 
criterion is used. We find a possibility that for m s < 200 MeV, the chiral broken vacuum 
phase might continuously connected with the CFL phase at high density without any phase 
transition. 

There are several future problems to be studied. 

Dirac mass function and chiral condensate : In order to study the phase boundary 
region where the chiral restoration occurs, we have to include the diagonal self-energy in 
our SD equation. Within the 4-Fermi contact interaction model, the competition between 
the chiral condensate and diquark condensate is investigated by several authors 33 )> 34 ). Also, 
the SD approach is examined in Ref. 20 ) for m s = oo. In particular, the ss condensate, or 
equivalently, the dynamical strange quark mass is shown to play an important role for the 
phase boundary between the 2SC and the CFL 10 ), which strongly stabilizes the 2SC phase 
with finite strange quark mass. The critical strange quark mass for unlocking transition may 
be significantly lowered by this effect 10 ^' 12 K 
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Wavefunction renormalization and chemical potential renormalization : There is an 
infrared singular factor in the fermion wavefunction renormalization : Z^ 1 (ko) = 1 + 
const. g 2 log(/z/<7o) at the one loop level 37 **. Furthermore, the chemical potential also suf- 
fers from a renormalization because if \x 7^ 0, Lorentz symmetry is explicitly broken down 
to 0(3). The full self-consistent treatment of both the gap, Z~ x and a chemical potential 
renormalization fi r is needed for more quantitative argument. 

Exact treatment of the strange quark mass : In this paper, the effect of the strange 
quark mass was examined within a simple kinematical criterion. It would be interesting to 
examine the full treatment at finite temperature in our improved SD approach 38 \ So far, 
the unlocking phase transition at finite temperature has not yet studied in a perfect way. 
Including the dynamical effect of the strange quark mass on the gaps, we can make sure 
whether the critical end point might show up at finite chemical potential or not, namely the 
end-point obtained in Ref. 12 -* is an artifact of an effective model with cut-off or not. 

Electric charge and color neutrality : Imposing the color and the electric charge neutrality 
has been shown to strongly disfavor the 2SC by the model independent analysis by Alford 
and Rajagopal 39 - 1 . Their argument based on the weak coupling expansion of thermodynamic 
potential in A ~ m%/f/, reveals that the 2SC state hardly realizes in the neutron stars. 
Several works 40 - ) ~~ 42 - 1 also show that the naive 2SC phase costs the large energy for enforcing 
the neutrality. It would be interesting to examine our improved SD approach to the neutral 
quark matter, and to determine the phase diagram under neutral condition. 

Another possibility of the ground state : As long as we restrict ourselves to the homo- 
geneous system, the unlocking transition is unique for a given chemical potential //, and 
thus the m s dependent unlocking critical line appears in the phase diagram. However, if we 
remove this restriction, then this line might split into two lines and another inhomogeneous 
state like the LOFF state 43) ' 44) might intrude between the CFL and the 2SC. For the 2SC 
phase with the chemical potential difference 8/j, — /j, u — fid, this phenomena has been inves- 
tigated in Ref. 45 - ) . Recently, it has been suggested that the interior gap state may overcome 
the LOFF state if the coupling between heavy quark and light quark exceeds some critical 
coupling 46 -* . It would be interesting to study these other possibilities than the pure BCS 
state in our SD approach with the momentum-dependent improved coupling. 
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